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We study heterostructures of singlet superconductors (SC) and strongly spin-polarized ferromag- 
nets (sFM) and show that a relative phase arises between the superconducting proximity amplitudes 
in the two ferromagnetic spin bands. We find a tunable pure spin supercurrent in a sFM contacted 
with only one SC electrode. We show that Josephson junctions are most effective for a spin polariza- 
tion P ~ 0.3, and that critical currents for positive and negative bias differ for a high transmission 
Josephson junction, due to a relative phase between single and double pair transmission. 
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Superconductor(SC) / ferromagnet(FM)-hybrid struc- 
tures have triggered considerable research activities in 
recent years 0, 0, 0, 0, 0, 0, 0, 0, 0, EE [H|. In Par- 
ticular, FM Josephson junctions arc promising spintron- 
ics devices as they allow for tuning the critical current 
via the electron spin. However, due to the competi- 
tion between the uniform spin alignment in the FM and 
spin-singlct pairing in the SC, singlet superconducting 
correlations decay in the FM on a much shorter length 
scale than in a normal metal [13]. Although this results 
in a rapidly decaying Josephson current for long junc- 
tions, the proximity effect leads to interesting physics 
in short and/or weakly polarized junctions, e.g., oscil- 
lations of the supercurrent as a function of the thick- 
ness of the interlayer that can give rise to 7r-junction 
behavior [13 , Il3j |. Recently however, in contradiction 
with these expectations, long-range supercurrents have 
been reported through strongly spin-polarized materials 
Q . Theoretical calculations have shown that for strongly 
polarized ferromagnets (sFM) spin scattering at SC/FM 
interfaces [1J] leads to a transformation of singlet corre- 
lations in the SC into triplet correlations Q (the 'triplet 
reservoirs' of Ref. 0]), that can carry a long-range su- 
percurrent through the sFM 0, 0, 0, 0, 0, [llj . 

So far, transport calculations in SC/FM hybrids have 
mostly concentrated on either fully polarized FMs, so- 
called half metals (HM), or on the opposite limit of 
weakly polarized systems. However, most FMs have an 
intermediate exchange splitting of the energy bands of 
the order of 0.2-0.8 times the Fermi energy E-p. For this 
intermediate range, one could naively expect a behavior 
similar to two shunted half metallic junctions. We will 
show, using a microscopic interface model, that this pic- 
ture is inadequate, and point out the crucial role played 
by the interfaces in coupling the sFM spin bands. 

In this work, we study Josephson junctions with a 
strongly polarized interlayer, and find fundamental dif- 
ferences compared to both half metallic and weakly po- 
larized interlayers. In particular, we see that although 
correlations between |- and |-electrons are suppressed 



due to the strong exchange field, spin-active interfaces 
generate interactions between long-range triplet super- 
currents in the two spin bands. We find that the long- 
range critical Josephson current varies non-monotonously 
with spin polarization P, showing a maximum around 
P = 0.3. Furthermore, additional phases arising from 



the interfaces [l4 |. specifically when the exchange split- 
ting is strong, lead to different current-phase relations 
for the spin-resolved currents I-\ and Iy through the junc- 
tion. We show how this gives rise to (i) a relative phase 
between single pair and "crossed" two-pair transmission; 
the latter process is illustrated in Fig. [T^, with equal 
numbers of pairs transferred in the spin ] and spin j 
band; (ii) different critical Josephson currents for oppo- 
site bias; (iii) equilibrium shifts in the current phase re- 
lation, in contrast to previous predictions @; and (iv) a 
tunable spin supercurrent in a FM brought into contact 
with a single SC electrode; we propose an experiment to 
measure this remarkable effect. 

Quasiclassical Green's functions (QCGF) [HI, [l6| are a 
powerful tool to describe hybrid structures of supercon- 
ductors and non-superconducting materials. Consider, 
e.g., the interface between a SC and a sFM shown in 




FIG. 1: (color online) (a) The coherent transfer of singlet 
pairs via a sFM (top) is not possible. However, the "crossed" 
pair transmissions process (bottom) is possible and leads to 
intriguing effects in high transmission junctions, (b) SC/sFM 
interface, showing the Fermi surfaces on either side (thick 
lines). Assuming momentum conservation parallel to the in- 
terface (fc||), a quasiparticle incident from the SC can either 
scatter into two (dotted arrows), or into only one (dashed 
arrows) spin band of the FM. 
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Fig. [TJa. For trajectories on the SC side, labeled 1, and 
characterized by Fermi momentum p-p\ and Fermi veloc- 
ity tTpi) the QCGF is obtained from the microscopic one, 
G, by integrating out the components oscillating on the 
Fermi wave length scale Xfi = fi/pFi- 9(pFi, R, £, t) = 
J d£, p T 3 G(p, R, e, t), where £ p = v F i(p-p F i). The QCGF, 
g, then varies as a function of the spatial coordinate R 
at a scale set by the superconducting coherence length 
£o = hv-Fi/2wkBT c , and obeys the Eilcnberger equation 



ihv F1 ■ V + [ef 3 — A — h, g] = 6, 



(1) 



with normalization condition g 2 = — 7r 2 1 Here, the 
hat denotes the 2x2 Nambu matrix structure in particle- 
hole space and t 3 is the third Pauli matrix; h includes all 
mean field and self energy terms governing the quasipar- 
ticle motion along QC trajectories aligned with v-px, and 
labeled by p F \] A is the SC order parameter. 

The exchange field Jfm in a sFM is comparable to 
the Fermi energy. As opposed to the weak polariza- 
tion limit (Jfm << E F ) this cannot be described by a 
term — Jfm ■ <? (with a the vector of Pauli spin matri- 
ces) in h of Eq. |T]), because the QC approximation in 
this case neglects terms of order Jp M /E F compared to 
A. In most SCs, this is not justified for Jfm > 0.1£f- 
However, for sufficiently large Jfm 3> y/EpA the coher- 
ent coupling of the spin bands in the FM can be disre- 
garded. Consequently, we define an independent QCGF 
for each spin band r\ € {2, 3} in Fig. Qi>: g(p Frn R, e, t) = 
J d£ pv T 3 G(p,R,s,t), where £ p) , = v F>1 (p-p Fv ). The ex- 
change field is incorporated by the different Fermi veloci- 
ties v Fri and momenta p Fr) in the two spin bands, and does 
not enter the equation of motion ([l} for the QCGFs. The 
g are Nambu matrices with diagonal (<?) and off-diagonal 
(/) components. These components are spin scalar, as 
opposed to the QCGF in the SC where they form a 2 x 2 
spin matrix as a result of spin coherence. Indeed, the 
spins of the pair wavefunction in the FM arc fixed cither 
to |tt) (band 2 in Fig-EEJa) or to |||) (band 3). 

The interface enters the QC theory in the form of ef- 
fective boundary conditions 
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I connecting the 
incident and outgoing QCGFs for the three Fermi surface 
sheets r\ e {1,2,3}. The boundary conditions are sub- 
ject to kinetic restrictions [2(J, as illustrated in Fig. [Tb. 
Note that, for a sFM, all singlet correlations are de- 
stroyed within the interface region (they decay on the 
short length scale X j = h/(p F2 -p F3 ) < ^«F2,3/A = £or, 
j2l|). The boundary conditions are formulated in terms 
of the normal-state scattering matrix (S-matrix) of the 
interface [l9[ , which for three bands has the general form 



S = ip 
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We obtain the reflection and transmission coefficients 
from a microscopic calculation. Wc consider an interface 



formed by a thin (w X F ) insulating FM layer of thick- 
ness d between the SC and bulk sFM (yellow areas in 
Figs.[T|D and[2k), characterized by an interface potential 
Vj — 3\ ■ 3. The orientation of the exchange field J\ in 
the interface layer is determined by angles a and tp, with 
a the angle between J\ and the exchange field Jfm of 
the bulk sFM (see Fig.[2}3). The S-matrix connecting in- 
and outgoing amplitudes in the bulk SC and sFM is then 
obtained by a wave-matching technique, where the am- 
plitudes in the interface layer are eliminated. Doing so, 
we obtain in the tunneling limit an S-matrix of the form 
Ru = e iW2)° 3 f 12 = f 2 i = (t2e«»/ 2 ,f 2 e- Wa / 2 ) r 
and f 13 = f 31 = (t' 3 e M3 / 2 ,t 3 e- M3 / 2 V. The spin mixing 
fJ-angles in these expressions [H [Til [Til (also called spin- 
dependent interfacial phase shifts [23j), and all remaining 
S-matrix parameters are obtained from a microscopic cal- 
culation as outlined above. As such, they depend on d, 
V\, a , tp, and the Fermi-momenta of the three bands 
(we assume |Ji| = | Jfm|)- The dependence on the an- 
gle tp is made explicit in Eq. while the dependence 
on the angle a is implicit in the r and t parameters via 
t' 23 oc sin(a/2), £2,3 oc cos(a/2) and r 23 , r 32 oc sina. In 
the following we use these tunneling-limit expressions to 
gain insight into the physics of the problem. The results 
shown in the figures however, are obtained by a full nu- 
merical calculation. For definiteness, we present results 
for parabolic electron bands with equal effective masses. 

Applying these boundary conditions to a Josephson 
junction depicted in Fig. [2^, and assuming bulk solutions 
for the QCGFs incoming from the SC electrodes, we ar- 
rive at the following system of linear equations for the 
/-functions in the tunneling limit (labels k,j e {L,R} 
with j k denote the left / right interfaces) : 



(3) 



Here, the factors (3 V = e - 2 \ 6 ™\ L / v ±n ; where L Xj is 
the junction length, e n = (2n + l)irk-QT the Matsubara 
frequency, v± v the Fermi velocity component along the 
interface normal, and 77 € {2, 3} the band index, arise 
from the decay of the / functions in the sFM layer. As 
depicted in Fig. [2^, coupling between the sFM spin bands 
is provided by the quantity (for our model r2 3 r 32 is real) 
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while the inhomogencity in Eq. (|3|), [j4i i; ]j, can be inter- 
preted as a pair transmission amplitude from the SC into 
spin band 77 of the sFM through the interface j. It reads 



[A 



sgn(£„) 



1-S 2 
B ri = r+/n 
sin?^^ ± s'm(i} v — 1?) 



{B v + C„)v;Ae l(X± ^ 



5 = A-sin(i?/2)/0^ 



(5) 

(6) 
(7) 



3 




a l \ : \ :a r 

i t i \v 



0.2 0.4 0.6 0.8 1 

k,/k 




(a) 



> 



A(p=0 






■ ■■■I, 

1 


- A(p=7t/4 









0.5 4 1, 1.5 
Ax/7C 




1 d o Acp/n i 



FIG. 3: (color online) (a) Spin resolved and total CPR for 
Atp — and 7r/4. (b) Coefficients I cp and Jo CT of Eq. © 
vs. Ax- (c) Critical current in positive(/ c +)/negative(7 c _) 
bias direction vs. Atp. (d) The equilibrium phase difference 
Axcq vs. Atp varies from n to 0. In all plots T — 0.2 T c , 
d — 0.25Xfi, P — 0.21, other parameters as in Fig. [5] 



FIG. 2: (color online) (a) Josephson-j unction with spin-active 
SC/sFM interfaces formed by magnetized layers (yellow), (b) 
Orientation of the interface magnetization described by spher- 
ical angles a and tp. (c) The quantities in Eq. ([7]) vs. fen for 
two polarizations P, and vs. P for perpendicular impact (in- 
set), (d) Critical current I c vs. temperature T for various po- 
larizations P of the sFM-layer. (e) 7 c 7? n -product and normal 
state resistance R n A as function of P for T = 0.5T C , d = Xyi, 
and (Vi - Ji)/Ef = 10 -4 (dotted), 0.2 (full), 0.5 (dashed). 
R^A in units of (e 2 iVpi«pi) _1 , Nfi being the normal state SC 
DOS. A = 1.76 meV. In all plots: «l — an. = 7r/2, tp-L = tpR, 
L = £o, d = 5Xfi, Vi — Ji = 0.5 Ef, pf2 = 1.18 Pfi, unless 
stated otherwise. P is tuned by j>F3. 



where f2 n = + A 2 , x is the order parameter phase 
of the corresponding SC, and the +(— ) sign in x^f cor ~ 
responds to 77 = 2(3). Note that t n t' n oc sin(a), implying 
that the generation of triplet correlations relics on a 7^ 0. 

In Fig. [2fc we show r^, Eq. for the majority spin 
band (?y = 2) as a function of fen. For large enough fen, 

T2 vanishes in contrast to '. This region of fen -values 
allows for transmission into only a single spin band of 
the FM (see Fig. [TJd). With increasing spin polarization 
P = (PF2 - PF3)/{PF2 + PF3) it extends over a larger 
range of fc|| -values, eventually spanning the entire Fermi- 
surface for a half-metal. At the same time the maximal 
value of t 2" decreases to zero, as demonstrated in the 
inset in Fig. [2^, where the parameters are shown for 
normal impact as function of P. The enter the B v and 
terms in Eq. ([6]), which exhibit different temperature 
(T) dependencies due to the additional \s n \ term in C n 
[7| ■ This interplay leads to an intriguing change in the In- 
dependence of the Josephson current, plotted in Fig. [2ji. 
For high P, a non-monotonic behavior is observed similar 
to that for a HM 0, |Tl| , due to the dominant C2 term, 
whereas for smaller P the term arising from P>2 leads to 
a monotonic decay with increasing T. As a result, the 
bump in I C (T) disappears with decreasing polarization. 



In Fig. [5^ we plot the I c R n product as a function of P 
(left scale). The variation of the normal state resistance 
R n with P (right scale) cannot account for the variation 
of I c . The critical current is suppressed for small P due 
to small spin mixing angles (see Fig. [5£), and for high 
P due to reduction of conductivity in the minority spin 
band. We thus predict a maximum critical current in a 
sFM junction for intermediate P ~ 0.3. We caution that 
in the hatched regions in Fig. [2^, there arc additional 
processes, not included in our model; e.g., for small P 
spin coherence leads to singlet amplitudes in the FM. 

We now discuss intriguing effects associated with the 
angles ^l,r [sec Eqs. (j4][5]) and Fig.|2b]- In Fig. [3k w e plot 
the spin-resolved current-phase relation (CPR) 24| for a 
high transparency junction (d = 0.25 Xpi) as a function 
of Ax = Xr ~ Xl for two values of Atp = tpn — ipi, [26|. 
Clearly, there is a non-trivial modification of the CPR in 
the presence of Aip. We find that the CPR can be well 
described by the leading Fourier terms in Atp, 



I cp - loa ■ sin ( Ax CT + oAtp 



(8) 



where a = +(— )1 for spin t(D- Here, 7o<r (shown 
in Fig. [3b) and Axo- are renormalized due to multiple 
transmission processes. The first term in Eq. ([8|) de- 
scribes a special type of multiple transmission process, 
which wc call "crossed pair" (cp) transmission, shown in 
Fig. [Tk. It is a result of singlet-triplet mixing and triplet 
rotation induced by the interfaces. Here two singlet 
Cooper pairs are effectively recombined coherently into 
two triplet pairs that propagate in different spin bands. 
Similar processes recombining a higher (but even) num- 
ber of pairs will also contribute. The phase associated 
with these processes comes from [^4i2^4i3]l[-^12-^13]r fac- 
tors with Ai v from Eq. ([5|), and is given by multiples of 
(Ax + Atp) + (Ax - Atp) = 2Ax- Consequently, 7 cp 
is independent of Atp and 7r-periodic in Ax, as shown 
in Fig. (foh li nc ) ■ It is also obvious that I cp is spin 
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FIG. 4: (color online) (a) Setup with only one SC electrode, 
(b) Spin-supercurrent I s vs. P. for various (Vi — J\)/E-p = 
10" 4 (dotted), 0.2(full), 0.5(dashed). R n A refers to the nor- 
mal state resistance of the SC/FM interface, d = Xfi, other 
parameters as in Fig. [2] 



symmetric, i.e., it carries a charge current, but no spin 
current. We find that transfer processes with even num- 
ber of pairs, but non-zero total spin, are in contrast to 
the cp transmission strongly suppressed. Contributions 
to the second term in Eq. (jHJ) come from processes that 
transmit one more Cooper pair in one of the spin bands 
compared to the other, including single pair transmis- 
sion. It is therefore spin-dependent in magnitude (see 
Fig. [3Jd) and shows Aip phase shifts with opposite signs 
for opposite spins. The relative phase between the two 
terms in Eq. ([8]) leads to surprising measurable effects 
for finite Aip and intermediate P. First, we find a dif- 
ference in the positive (I c +) and negative (I c —) bias crit- 
ical charge currents, as shown in Fig. [3fc. This is also 
directly visible in Fig. [3^, where the maximum and min- 
imum current have a different absolute value. Second, 
we find a shift of the equilibrium phase A% C q for the 
charge current, as shown in Fig. [3ji (the jump as func- 
tion of Aip is associated with multiple local free energy 
minima). We note that in the tunneling limit, Eq. ([8|) re- 
duces to I a w — Iqo- ■ sin(Ax + cr Aip) , and the equilibrium 
phase shift is present as long as Jot 7^ ^04. • 

Another remarkable consequence of a non-zero Ap is 
observed for a setup shown in Fig. [4^, when an sFM is 
coupled via a spin-active interlaycr to a single SC on the 
left, and is terminated by a magnetic surface on the right. 
All quasiparticles are reflected at the surface, leading to 
a zero charge current. However, not all quasiparticles are 
reflected back into their original spin band since spin-flip 
reflections [^23 in Eq. Q] mediate interactions between 
the two spin bands, and, remarkably, a pure spin su- 
percurrent remains. In this case, both terms in Eq. (|8|) 
vanish as they are related to direct transmission. In- 
stead, the leading term for the spin supercurrent is of 
second order in Ap, I oc sin(2Ay>), resulting from the 
phases picked up when a triplet Cooper pair reflects at 
the right interface [25| . The maximal spin-current, de- 
fined as I s = max.A.<pI(A(p), is plotted in Fig. [4}d as a 
function of spin polarization. Note that it vanishes both 
for P — > and P — ► 1, since it requires the presence of 
two bands, and is maximum for intermediate P. 

This pure spin current can be tuned by an external 



microwave field that couples to the magnetization of the 
right surface in Fig. [4ji, and thus leads to a time depen- 
dent Ap(t). A high degree of control can be achieved 
by manufacturing a surface layer using a different mag- 
netic material, preferably magnetized perpendicular to 
the bulk FM, thus optimizing external tunability. As 
Ap(t) acts as a time dependent superconducting phase, 
we predict in addition to a spin accumulation in the FM 
a measurable ac spin supercurrent, analogously to the ac 
charge Josephson current in a voltage biased junction. 

In summary, we have presented a study of heterostruc- 
tures between singlet superconductors and strongly spin- 
polarized ferromagnets. We have found that the Joseph- 
son effect markedly differs from that for a fully polar- 
ized material or for a ferromagnet with a weak spin band 
splitting. We discussed the importance of the phase-shift 
between single pair and "crossed" two-pair transfer pro- 
cesses that leads to measurable anomalous junction be- 
havior. We have also found that a pure spin supercurrent 
is induced in a strongly polarized ferromagnet coupled to 
one singlet superconducting electrode, and have proposed 
a way of measuring this effect. 

We thank T. Lofwander for stimulating discussions. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[«] 
[9] 
[10] 

[11] 
[12] 

[13] 

[14] 
[15] 

[16] 

[17] 

[18] 
[19] 



[20] 

[21] 
[22] 
[23] 



V. 
Y. 
M. 
V. 



F. S. Bergeret et al, Rev. Mod. Phys. 77, 1321 (2005). 
A.I. Buzdin et al, Rev. Mod. Phys. 77, 935 (2005). 
M. Eschrig et al. Phys. Rev. Lett. 90, 137003 (2003). 
A.F. Volkov et al, Phys. Rev. Lett. 90, 117006 (2003). 
J. Kopu et al, Phys. Rev. B 69, 094501 (2004). 

R.S. Keizer et al., Nature 439, 825 (2006). 
M. Eschrig et al, J. Low Temp. Phys. 147, 457 (2007). 
M. Houzet et al, Phys. Rev. B 76, 060504(R) (2007). 
Braude and Yu. Nazarov, Phys. Rev. Lett. 98 (2007). 
Asano et al, Phys. Rev. Lett. 98, 107002 (2007). 
Eschrig et al, Nature Phys. 4, 138 (2008). 
Buzdin et al, Pis'ma Zh. Eksp. Teor. Fiz. 35, 147 
(1982); Z. Radovic et al, Phys. Rev. B 68, 014501 (2003). 
V.V. Ryazanov et al, Phys. Rev. Lett. 86, 2427 (2001); 
T. Kontos et al, Phys. Rev. Lett. 89, 137007 (2002). 
T. Tokuyasu et al, Phys. Rev. B 38, 8823 (1988). 
A. Schmid and G. Schon, J. Low Temp. Phys. 20, 207 
(1975); J.W. Serene et al, Phys. Rep. 101, 221 (1983). 

G. Eilenberger, Z. Phys. 214, 195 (1968); A.I. Larkin and 
Y.N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 55, 2262 (1968). 
A. Shelankov, J. Low Temp. Phys. 60, 29 (1985); A. 
Shelankov and M. Ozana, Phys Rev. B 61, 7077 (2000). 
M. Eschrig, Phys. Rev. B 61, 9061 (2000). 

A. Millis, D. Rainer, and J.A. Sauls, Phys. Rev. B 38, 
4504 (1988); M. Fogelstrom, ibid. 62, 11812 (2000); E. 
Zhao et al, ibid. 70, 134510 (2004). 

I. Zutic and O.T. Vails, Phys. Rev. B 61, 1555 (2000); 

K. Halterman and O.T. Vails, ibid. 66, 224516 (2002); J. 

Linder and A. Sudb0, ibid. 75, 134509 (2007). 

The condition for treating the singlet components within 

2 2 
QC theory is Xj » CotjA'f,, » X ¥v . 

With a proper choice of the quantization axis in the SC, 
the reflection matrix Rn is spin-diagonal. 
A. Brataas et al, Phys. Rev. Lett. 84, 2481 (2000); D. 
Huertas- Hernando et al, ibid. 88, 047003 (2002); A. Cot- 
tet and W. Belzig, Phys. Rev. B 72, 180503 (2005). 



5 



[24] A. A. Golubov el ai, Rev. Mod. Phys. 76, 411 (2004). [26] Our final results only depend on Aip = tp R - <£l, in ac- 

[25] A similar effect, however involving triplet superconduc- cordance with global spin rotational invariance. 

tors with a FM interlayer has been reported in P. M. R. 

Brydon, D. Manske larXiv:0901.4096l 



